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ABSTRACT:  Graphene-based  devices  have  garnered  tremendous 
attention  due  to  the  unique  physical  properties  arising  from  this  purely 
two-dimensional  carbon  sheet  leading  to  tremendous  efficiency  in  the 
transport  of  thermal  carriers  (i.e.,  phonons).  However,  it  is  necessary 
for  this  two-dimensional  material  to  be  able  to  efficiently  transport 
heat  into  the  surrounding  3D  device  architecture  in  order  to  fully 
capitalize  on  its  intrinsic  transport  capabilities.  Therefore,  the  thermal 
boundary  conductance  at  graphene  interfaces  is  a  critical  parameter  in 
the  realization  of  graphene  electronics  and  thermal  solutions.  In  this 
work,  we  examine  the  role  of  chemical  functionalization  on  the 
thermal  boundary  conductance  across  metal/graphene  interfaces. 

Specifically,  we  metalize  graphene  that  has  been  plasma  functionalized 
and  then  measure  the  thermal  boundary  conductance  at  Al/ graphene/ 

Si02  contacts  with  time  domain  thermoreflectance.  The  addition  of  adsorbates  to  the  graphene  surfaces  are  shown  to  influence 
the  cross  plane  thermal  conductance;  this  behavior  is  attributed  to  changes  in  the  bonding  between  the  metal  and  the  graphene, 
as  both  the  phonon  flux  and  the  vibrational  mismatch  between  the  materials  are  each  subject  to  the  interfacial  bond  strength. 
These  results  demonstrate  plasma-based  functionalization  of  graphene  surfaces  is  a  viable  approach  to  manipulate  the  thermal 
boundary  conductance. 
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Graphene-based  devices  have  garnered  tremendous  atten¬ 
tion  due  to  the  unique  physical  properties  arising  from 
this  purely  two-dimensional  carbon  sheet.1,2  These  unique 
properties  allow  for  tremendous  efficiency  in  the  transport  of 
not  only  charge  carriers  but  thermal  carriers  (i.e.,  phonons)  as 
well.3”  For  these  reasons,  graphene  is  being  pursued  both  for  a 
variety  of  electronic  applications  and  as  an  enabler  of  next- 
generation  thermal  solutions.  As  such,  it  is  necessary  for  this 
two-dimensional  material  to  be  able  to  efficiently  transport  heat 
into  the  surrounding  3D  device  architecture  in  order  to  fully 
capitalize  on  its  intrinsic  transport  capabilities.  Therefore, 
efficient  heat  flow  across  solid-graphene  interfaces  as  described 
by  the  thermal  boundary  conductance,8  or  Kapitza  con¬ 
ductance,9  hx,  is  vital  for  the  realization  of  graphene  electronics 
and  thermal  solutions. 

To  this  end,  previous  groups  have  measured  hK  at  Si02/ 
graphene  (ref  10)  and  Au/graphene/SiOz  interfaces  (ref  4). 

The  Si02/graphene  interface  exhibited  a  thermal  boundary 
conductance  at  room  temperature  of  ~90  MW  m-2  IC1  (ref 
ll).  In  contrast,  the  Au-graphene  interface  demonstrated  a 
thermal  boundary  conductance  at  room  temperature  that  was 
nearly  a  factor  of  2  lower  (refs  4  and  12).  This  discrepancy  can 
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be  explained  utilizing  traditional  phonon  scattering  models 
(e.g.,  diffuse  mismatch  model,  or  DMM  )  in  which  the  lower 
conductance  at  the  Au  boundary  is  attributed  to  this  material 
having  a  lower  energy  cutoff  of  available  modes  available  to 
transport  energy  as  compared  to  both  the  graphene  and  Si02. 
While  this  theory  is  supported  in  part  by  investigations  of 
metal/graphite  interfaces,  chemical  bonding  between  the  layers 
has  been  shown  to  play  a  prominent  role  in  phonon 
transmission  and  thus  hK  as  well.  4  In  spite  of  this  fact,  the 
interplay  between  chemical  bonding,  the  allowable  phonon 
modes,  and  the  subsequent  effect  on  interfacial  thermal 
transport  remains  unclear.  For  graphene  boundaries,  mean¬ 
while,  the  material’s  inherent  surface  sensitivity  will  only 
enhance  this  interplay. 

In  response,  we  investigate  the  role  of  interfacial  bonding  at 
metal/single-layer-graphene  (SLG)  interfaces  by  introducing 
chemical  adsorbates  on  the  SLG  surfaces  in  order  to  increase 
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the  density  of  covalent  bonds  bridging  the  metal  and  the  SLG. 
In  doing  so,  we  demonstrate  plasma-based  functionalization  of 
graphene  surfaces  as  a  means  to  manipulate  the  thermal 
boundary  conductance.  Specifically,  we  metalize  plasma- 
functionalized  graphene  and  then  measure  hK  at  Al/SLG/ 
Si02  contacts  with  time  domain  thermoreflectance 
(TDTR). 15,16  From  these  measurements,  adsorbates  on  the 
SLG  surfaces  are  shown  to  influence  the  cross  plane  thermal 
conductance.  This  is  attributed  to  changes  in  the  bonding 
between  the  metal  and  the  SLG.  Additionally,  the  Al/SLG/ 
Si02  thermal  boundary  conductance  is  shown  to  consist  of  two 
individual  conductances,  Al/SLG  and  SLG/Si02,  in  line  with 
the  conclusions  from  ref  4.  Using  our  measured  values  of  hK  in 
conjunction  with  previous  measurements  between  metals  and 
carbon-based  materials,4’14  the  thermal  boundary  conductance 
is  found  to  be  most  influenced  by  the  interfacial  chemical 
bonding  as  both  the  phonon  flux  and  the  vibrational  mismatch 
between  the  materials  are  each  subject  to  the  interfacial  bond 
strength. 

The  graphene  films  were  grown  by  chemical  vapor 
deposition  on  Cu  substrates  and  transferred  to  Si02/Si 
substrates  using  the  conventional  wet  chemical  approach.17 
The  graphene  films  were  then  functionalized  with  either  oxygen 
or  hydrogen  using  electron  beam  generated  plasmas18 
produced  in  gas  mixtures  containing  02  and  H2,  respectively. 
Electron  beam  generated  plasmas  are  particularly  attractive  in 
that  they  offer  high  plasma  densities  (lO10— 1011  cm-3)  along 
with  low  electron  temperatures  (<1  eV),  which  provides  large 
fluxes  of  low  energy  (<5  eV)  ions11  thereby  eliminating  ion- 
induced  etching  or  sputtering  of  the  basal  plane.  For  this  work, 
the  process  parameters  include  a  pressure  of  50  mTorr  Ar 
diluted  with  either  02  or  H2,  where  the  reactive  gas  was  5%  of 
the  total  flow  rate.  The  plasma  was  operated  at  a  duty  factor  of 
10%  (2  ms  pulse  width  and  20  ms  period)  and  the  total 
treatment  time  was  60  s  (plasma  exposure  time  was  6  s.).  The 
functionalization  process  leads  to  an  inclusion  of  chemical 
moieties  on  the  graphene  surface,  which  was  verified  by  X-ray 
photoelectron  spectroscopy  (XPS).  The  total  oxygen  present 
on  the  surface  after  treatment  in  the  oxygen  containing  plasma 
is  ~25  at.%.  A  careful  inspection  of  the  Cls  peak  reveals  the 
presence  of  hydroxyl  and  carboxyl  groups  as  indicated  by  the 
well-pronounced  peaks  around  286.6  and  289  eV  (Figure  la). 
Although  the  presence  of  hydrogen  cannot  be  quantified  by 
XPS,  the  broadening  of  the  main  peak  (located  at  384.5  eV) 
indicates  the  increase  of  sp3  hybridized  carbon  atoms,  which  is 
mainly  attributed  to  the  incorporation  of  hydrogen.  From  the 
fitting  procedure,  the  contribution  from  the  sp3  peak  (situated 
at  +0.65  eV  from  main  peak)  increased  from  22%  for  graphene 
to  ~33%  for  FI-treated  graphene  (Figure  la). 

To  structurally  characterize  the  graphene  films,  we  perform 
Raman  spectroscopy  on  the  functionalized  samples  along  with  a 
reference  film.  In  each  case,  a  Raman  spectrum  is  acquired 
every  333  nm  across  a  15  X  15  ^m  area  using  532  nm  laser  light 
resulting  in  a  data  set  of  2025  spectra  for  each  sample.  Laser 
powers  were  not  found  to  damage  the  sample.  Using  the  fwhm 
of  2D  mode  at  ~2678  cm-1,  each  of  the  films  are  found  to  be 
composed  almost  completely  of  SLG.20  As  very  little  difference 
in  the  character  of  the  graphene  is  found  across  the  sample, 
each  of  the  SLG  spectra  are  averaged  and  shown  in  Figure  lb. 

Before  functionalization,  the  graphene  films  contain  very 
little  disorder  as  the  ratio  of  the  D-band  (~1350  cm-1)  intensity 
to  that  of  the  G-band  (~1580  cm-1)  is  I(D)/I(G)  «  0.2. 
Significant  levels  of  disorder  become  apparent,  however,  upon 
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Figure  1.  (a)  High  resolution  Cls  spectra  of  untreated  graphene  and 
plasma  treated  graphene.  The  inset  shows  the  oxygen  treated  sample 
along  with  peak  fitting  results,  (b)  Averaged  Raman  spectra  for  each  of 
the  films  analyzed  in  this  study. 

exposure  to  either  the  oxygen  or  hydrogen  plasma.  Such 
disorder  is  readily  apparent  in  each  film  through  the  increased 
intensities  of  several  spectral  features  (e.g.,  D,  D',  and  D  +  D' 
modes)  that  are  directly  indicative  of  an  alteration  from  an  ideal 
graphene  lattice.  While  both  functionalized  films  do  exhibit 
disorder,  the  extent  of  this  disorder  is  substantially  different.  In 
the  hydrogen  exposed  film  for  example,  the  2D  peak  and  the  G- 
peak  have  much  greater  intensities  than  that  of  the  D  +  D'  peak 
and  the  D'  peak,  respectively.  In  contrast,  the  G  and  2D  features 
have  similar  intensities  in  comparison  to  the  D'  and  D  +  D' 
modes  for  the  film  exposed  to  an  oxygen  plasma.  In  total,  these 
results  indicate  increased  disorder  with  either  functionalization 
procedure  albeit  with  a  heightened  disruption  from  exposure  to 
the  oxygen  plasma  relative  to  that  of  hydrogen. 

The  origin  of  this  increased  defect  response  can  evolve  from 
both  plasma  induced  damage/amorphization  of  the  graphene 
or  changes  in  the  carbon— carbon  bonding  structure  that  occur 
with  covalent  bonding  of  the  oxygen/hydrogen  onto  the 
graphene  lattice.  Due  to  the  low  energy  of  the  plasma 
functionalized  procedure,  no  etching  or  changes  in  the  lattice 
are  expected  to  occur.18'”1  Therefore,  we  attribute  the  increased 
defect  response  of  the  Raman  signal  to  result  from  the  presence 
of  the  functionalized  species  (i.e.,  hydrogen/oxygen)  rather 
than  the  inducement  of  any  damage  into  the  graphene  lattice. 
Finally,  it  is  of  note  that  these  changes  do  not  alter  the  strain 
state2  or  carrier  concentration23  between  the  films  as  each  are 
found  to  have  peak  positions  of  the  G  and  2D  bands  that  are 
within  the  spectrometer  resolution  of  ±1  cm-1.  Thus,  any 
induced  changes  in  the  heat  transport  are  a  consequence  of  the 
bonding/defect  states  of  the  film  and  not  due  to  strain  or 
carrier  effects. 

Prior  to  TDTR  measurements,  we  metalize  each  sample  with 
90  nm  of  electron  beam  evaporated  A1  at  base  pressures  no 
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greater  than  2.3  X  10-7  Torr.  We  then  measure  the  thermal 
boundary  conductance  across  the  Al/ SLG/ Si02  interfaces  with 
TDTR24  TDTR  and  appropriate  analyses  accounting  for  pulse 
accumulation  when  using  a  Ti:  Sapphire  oscillator  have  been 
discussed  by  several  groups  previously.  '  ’  Our  specific 
experimental  setup  is  described  in  detail  elsewhere.16  In  short, 
TDTR  is  a  pump— probe  technique  in  which  our  laser  pulses 
emanate  from  a  Spectra  Physics  Mai  Tai  with  90  fs  pulses  at  80 
MHz  repetition  rate.  We  delay  the  time  in  which  the  probe 
pulse  reaches  the  sample  with  respect  to  the  pump  pulse  with  a 
mechanical  delay  stage  that  gives  ~4  ns  of  probe  delay.  We 
modulate  the  pump  path  at  11  MHz  and  monitor  the  ratio  of 
the  in-phase  to  out-of-phase  signal  of  the  probe  beam  from  a 
lock-in  amplifier  (— !/„/!/,  llt).  Our  pump  and  probe  spots  are 
focused  to  15  /fm  radii  at  the  sample  surface;  at  these  sizes 
when  modulating  at  1 1  MHz,  we  should  be  negligibly  sensitive 
to  any  in-plane  transport  in  the  graphene  sheet4,16  decreasing 
the  uncertainty  associated  with  determining  hK.  We  test  each 
sample  in  a  liquid  nitrogen  cooled  cryostat  and  measure  hK 
from  100  to  400  K.  We  perform  several  TDTR  scans  at 
different  locations  of  two  different  samples  of  each  sample  type 
(7—10  scans  total  for  each  type  of  graphene  sample).  Typical 
uncertainties  due  to  the  different  locations  on  the  samples  or 
from  sample  to  sample  were  less  than  10%  (calculated  from  the 
standard  deviation  among  all  the  TDTR  scans  on  a  given 
sample  type),  which  is  smaller  than  the  variation  in  metal/ 
graphene/Si02  hK  expected  due  to  a  transition  from  SLG  to  n- 
layer  graphene  where  n  >  1  (~25— 35%),  as  measured  via 
TDTR4  This  further  indicates  that  our  samples  were  primarily 
SLG,  confirming  our  Raman  analysis.  We  fit  the  TDTR  data 
with  a  model  that  accounts  for  pulse  accumulation  in  a  three- 
layer  system  (90  nm  Al,  275  nm  Si02,  and  Si).  We  adjust  hK 
between  the  Al  and  the  Si02  to  determine  the  Al/SLG/Si02 
thermal  boundary  conductance. 

Figure  2  shows  the  thermal  boundary  conductance  of  the 
various  graphene  samples  as  a  function  of  temperature,  T. 
There  is  a  clear  increase  in  the  overall  hK  of  the  Al/SLG/Si02 
contact  for  the  O-functionalized  graphene,  yet  a  slight  decrease 
in  hK  for  the  H-functionalized  graphene  as  compared  to  the 
reference  sample.  At  room  temperature,  we  show  an  increase  in 
hK  of  ~25%  due  to  O-fiinctionalization  compared  to  the 
nonfunctionalized  sample.  As  evident  from  the  XPS  and  Raman 
analysis,  both  the  H-  and  O-functionalization  leads  to  increased 
sp3  carbon  bonds.  In  hydrogen,  the  C— H  sp3  bond  is 
unreactive  and  thereby  leaves  the  graphene  surface  inert.28-30 
However,  it  is  well-known  that  the  addition  of  oxygen 
functional  groups  on  a  surface  greatly  changes  the  surface 
energy,  leading  to  an  enhanced  adhesion  at  a  metal/oxide 
interface.  This  is  mainly  due  to  the  oxide’s  electronic  and 
structural  properties.  It  was  shown  for  example,  that  graphene 
oxide  has  a  higher  surface  energy  compared  to  graphene.31 
Moreover,  in  the  C— O  bond,  one  electron  from  the  oxygen 
molecule  is  shared  with  the  carbon,  which,  in  turn  leaves  an 
additional  electron  to  bond  with  the  metal  surface,  thereby 
increasing  the  bond  strength  between  the  Al  and  the  graphene. 
We  also  note  that  both  the  H-SLG  and  O-SLG  exhibit  disorder 
as  compared  to  the  nonfunctionalized  SLG  (c.f.  Figure  lb). 
The  effects  of  disorder  on  the  cross  plane  thermal  conductance 
of  SLG  are  relatively  minor  and  lead  to  a  slight  decrease,  as 
determined  by  comparing  the  measured  hK  at  the  Al/SLG/Si02 
to  the  Al/H-SLG/Si02.  Unlike  our  previous  results  in  which  we 
had  shown  reductions  in  hK  due  to  roughness  at  Al/Si 
interfaces,32  the  functionalization  process  in  this  work  does  not 
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Figure  2.  Measured  thermal  boundary  conductances  on  the  Al/SLG/ 
Si02  sample  (filled  squares),  the  hydrogen  functionalized  sample 
(filled  circles),  and  the  oxygen  functionalized  sample  (filled 
diamonds).  The  hydrogen  functionalization  process,  which  introduces 
disorder  on  the  SLG,  also  leaves  the  SLG  surface  chemically  inert  and 
does  not  leave  any  additional  bonding  mechanism  for  the  Al  to  the 
SLG.  Functionalizing  the  SLG  with  oxygen  leaves  the  graphene  surface 
reactive  leading  to  increased  covalent  bonds  linking  the  Al  to  the  SLG, 
which  results  in  a  higher  phonon  transmission  and  increase  in  /zK.  For 
comparison,  we  also  show  hK  across  SLG/Si02  (open  circles,  ref  10), 
Au/Ti/SLG/Si02  (open  triangles,  ref  4),  and  Al/graphite  (open 
squares,  ref  14)  interfaces.  Using  these  two  measured  conductances 
with  eq  1,  we  can  predict  the  Al/SLG/Si02  interface  assuming  that  the 
Al/graphite  hK  is  similar  to  Al/SLG  and  that  the  Al/SLG  and  SLG/ 
Si02  conductances  can  be  separated.  The  resulting  conductance,  which 
is  depicted  by  the  solid  line,  is  in  very  good  agreement  with  our  Al/ 
SLG/Si02  measurements,  indicating  that  the  Al/SLG/Si02  thermal 
boundary  conductance  can  be  described  by  two  separate  conductances, 
Al/SLG  and  SLG/SiQ2. 


S  LG/SiC) 


o  o  o  o  o 


Al/O-SLG/SiO, 


Al/graphite 
\  Al/SLG/SiO, 


Au/Ti/SLG/SiO. 


AI/H-SLG/SiO, 


lead  to  a  significant  increase  in  the  surface  roughness  of  the 
graphene  films.  Therefore,  we  attribute  the  increase  in  hK  to  the 
enhanced  bonding  in  the  Al/0-SLG/Si02. 

For  comparison,  we  also  show  hK  measured  at  a  SLG/Si02 
interface,10  a  Au/Ti/SLG/Si02  interface  (the  Ti  is  a  2  nm 
wetting  layer),4  and  a  Al/graphite  interface.14  Using  these 
values  for  hK  at  the  Al/graphite  and  SLG/Si02  interface,  we  can 
compare  our  measured  conductance  across  the  Al/SLG/Si02 
interface  to  an  empirical  prediction  based  on  the  assumption 
that  the  overall  interface  conductance  can  be  separated  into  two 
conductances  in  series  on  either  side  of  the  graphene,  given  by 


h 


K, metal/ SI.G/SiO, 

1 


-1 


+ 


^  ^K, metal/ SLG  ^K,SLG/Si02 


(1) 


This  relation  assumes  that  the  cross  plane  thermal  resistance  of 
the  SLG  sheet  is  much  less  than  the  resistance  at  each  of  the 
interfaces,  which  is  a  valid  assumption  that  has  been  addressed 
in  detail  previously.4  eq  1  is  shown  in  Figure  2  as  the  solid  line, 
and  shows  relatively  good  agreement  with  our  measured  hK 
across  the  Al/SLG/Si02  interface. 

To  directly  compare  the  interface  conductance  between 
metals  and  SLG,  we  use  eq  1  to  estimate  hK  across  the  metal/ 
SLG  interface  by  assuming  that  the  SLG/Si02  thermal 
boundary  conductance  is  the  same  as  that  measured  by  Chen 
et  al.10  and  solving  for  li^metal/SLG  f°r  the  Al/SLG  and  Al/O- 
SLG  interfaces  in  this  work.  These  values  for  hK  at  the  Al/ SLG 
and  Al/O-SLG  interfaces  are  shown  in  Figure  3  along  with  hK 
Al/O-diamond  (ref  33)  Al/diamond  interfaces  (ref  34).  We  see 
a  factor  of  2  increase  in  hK  across  the  Al/O-SLG  interface  as 
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Figure  3.  Thermal  boundary  conductance  across  our  Al/SLG  (filled 
squares)  and  Al/ O-SLG  (filled  diamonds)  interfaces  derived  from  our 
experimental  measurements  across  the  Al/SLG/Si02  and  Al/ O-SLG/ 
Si02  interface  along  with  eq  1  and  the  data  from  the  SLG/Si02 
interfaces.10  The  thermal  transport  across  the  Al/SLG  interface 
increases  by  a  factor  of  2  with  oxygen  functionalization  of  the 
graphene.  Our  derived  values  for  /iK  across  the  Al/ SLG  interface  are  in 
good  agreement  with  conductances  measured  across  Al/graphite 
(filled  circles;  ref  14)  and  Al/diamond  interfaces  (open  squares;  ref 
34).  Similarly  our  derived  values  for  Al/O-SLG  hK  is  in  good 
agreement  with  Al/O-diamond  measurements  (open  diamonds,  ref 
33).  We  model  the  thermal  conductance  across  the  Al/SLG  interface 
with  the  DMM  (eqs  2—5),  as  shown  by  the  solid  line.  We  adjust  the 
velocity  of  the  SLG  in  our  DMM  calculations  to  model  the  Al/ O-SLG 
thermal  boundary  conductance  (dashed  line),  and  find  that  the 
resultant  velocity  in  the  oxygenated  SLG  is  about  a  factor  of  2  higher 
than  the  nonfunctionalized  sample.  This  is  indicative  of  the  increase  in 
bond  strength  between  the  Al  and  the  SLG  via  the  oxygen  adsorbates. 


compared  with  the  Al/SLG  interface.  We  attribute  this 
observed  behavior  to  the  enhanced  reactivity  of  the  oxygenated 
surface,  which  is  consistent  with  the  increase  observed  across 
the  Al/O-diamond  interface  as  compared  to  the  Al/diamond 
interface.  Recall  that  the  presence  of  oxygen-containing 
functional  groups  increases  the  surface  energy  and  thus 
makes  graphene  more  reactive,  allowing  a  lower  energy 
bonding  state  with  the  evaporated  Al  film  as  compared  to  the 
untreated  SLG  which  is  chemically  inert.28  30  This  lower 
energy  bonding  state  promotes  Al— O  formation,  leading  to  a 
stronger  Al— O— SLG  junction  and  thereby  an  increase  in  hK. 
Still,  it  is  interesting  that  while  hydrogen-termination  can 
significantly  decrease  the  conductance  of  Al/diamond  inter¬ 
faces,33  it  has  minimal  effect  on  Al/graphene  interfaces  (see 
Figure  2). 

To  quantify  this,  we  turn  to  the  diffuse  mismatch  model 
(DMM).8  Assuming  crystallographic  isotropy  in  the  Al  film,  the 
thermal  boundary  conductance  from  the  Al  to  the  SLG  is  given 
by 

=  A  2  f  dk 

8n  .  Jh  dT  (2) 

where  h  is  the  reduced  Planck’s  constant,  a>  is  the  phonon 
angular  frequency,  k  is  the  phonon  wavevector,  f  is  the  Bose— 


Einstein  distribution  function,  T  is  the  temperature,  v  is  the 
velocity,  and  f  is  the  phonon  transmission  coefficient  from  side 
1  to  side  2  (which  ranges  from  0  to  l).  Here,  the  subscript  “1” 
refers  to  the  corresponding  property  of  Al,  “2”  refers  to  the 
corresponding  property  of  graphene,  and;  refers  to  the  phonon 
polarization.  To  calculate  the  transmission  coefficient,  we  use 
the  assumption  discussed  by  Duda  et  al.35  in  which  we  treat  the 
graphene  sheet  as  a  two-dimensional  solid.  Therefore,  perform¬ 
ing  detailed  balance  on  the  fluxes  in  the  Al  and  SLG,36  we 
obtain 


?1- 


‘h 

<h  +  *2 


(3) 


where  q  is  the  phonon  flux.  For  the  flux  in  the  Al,  we  can  make 
the  usual  approximation  of  phonon  flux,37  which  is  the  basis  of 
the  derivation  of  eq  2, 38  namely 


i  =  A  X  /  H(fc)fc2Mfc)^2 dk 

8n  . 


(4) 


However,  given  the  two-dimensional  nature  of  the  SLG,  the 
phonon  flux  in  this  material  is  given  by35 

q2  =  -L  2  [  fcra,'(fc)fc>2,/fc)  dk 

Situ  .  Jk2  (5) 

where  a  is  the  interlayer  spacing  of  graphite  (i.e.,  a  =  3.35  A).39 
For  calculations  of  eqs  2—5,  we  fit  polynomials  to  the  phonon 
dispersions  of  Al  in  the  T  — >X  direction4  and  of  graphene  in  the 
r  — >  K  direction41  to  give  greater  accuracy  in  the  DMM 
calculations  as  compared  to  the  Debye  approximation.42  We 
note,  however,  that  the  modes  in  the  T  — >■  K  in  graphene  are 
extremely  Debye-like  (i.e.,  non  dispersive)  in  the  frequency 
regime  that  is  elastically  accessible  to  the  Al  (frequencies  below 
10  THz).  We  ignore  the  ZA  mode  in  the  graphene  dispersion 
since  the  ZA  mode  is  heavily  suppressed  in  encased  and 
supported  SLG.7,4j  We  note  that  our  calculations  in  eq  3  are 
balancing  the  total  flux,  and  not  the  flux  at  any  given  frequency. 
We  assume  elastic  scattering  in  this  calculation.  Therefore,  we 
restrict  the  wavevector  in  each  mode  of  the  SLG  and  perform 
the  integration  only  to  the  wavevector  of  the  corresponding 
mode.  As  a  result,  no  frequencies  above  the  cutoff  frequency  of 
Al  can  participate  in  hK.  Also,  given  that  our  measurements 
represent  the  cross  plane  conductance  of  Al/SLG,  the  cross 
plane  velocity  of  SLG  is  meaningless.  Therefore,  we  define  Vy 
as  a  fitting  parameter  that  we  assume  is  constant  with 
wavevector.  This  is  a  similar  approach  as  performed  in  the 
analysis  of  Koh  et  al.4  in  which  they  adjust  q2  to  fit  the  DMM  to 
their  Au/Ti/SLG  data  and  Schmidt  et  al.  in  which  they  adjust 
the  phonon  transmission  to  fit  the  DMM  to  their  metal/ 
graphite  data.14  We  note  that  this  approach  of  adjusting  only  Vy 
gives  us  more  direct  insight  into  how  the  bonding  at  the  Al/ 
SLG  interface  changes  due  to  the  functionalization  since  we  are 
not  adjusting  any  aspect  of  the  graphene  dispersion,  only  the 
transport  velocity.  The  fits  of  DMM  are  shown  in  Figure  3.  The 
solid  line  is  the  DMM  calculation  assuming  that  the  graphene 
cross  plane  velocities  are  v2L  =  2,455  m/s  and  v2T  =  1,480  m/s 
for  the  longitudinal  and  transverse  modes  of  the  SLG, 
respectively.  The  dashed  line  is  the  DMM  calculations 
assuming  v2jL  =  4,687  m/s  and  v2T  =  2,825  m/s  for  the 
longitudinal  and  transverse  modes  of  the  O-SLG,  respectively. 
The  velocities  that  result  in  the  best  fit  of  the  DMM  to  the  SLG 
data  are  in  good  agreement  with  the  cross  plane  velocity  of  bulk 
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graphite, 37  which  could  be  indicative  of  similar  bonding 
between  the  A1  and  SLG  as  the  van  der  Waals  bonds  cross 
plane  in  graphite.  The  velocities  resulting  in  the  best  fit  in  the 
Al/O-SLG  data  are  nearly  a  factor  of  2  higher,  which  we 
attribute  to  the  increased  bond  strength  between  the  A1  and 
SLG  due  to  the  presence  of  oxygen  leading  to  covalent  bonding 
with  the  A1  and  O-SLG.  It  is  important  to  note  that  we  use  the 
velocity  parameter  and  fitting  only  to  lend  insight  into  the 
nature  of  the  bonding  at  the  two  different  interfaces.  Since  the 
velocity  is  proportional  to  the  square  root  of  the  spring 
constant,  the  increase  by  a  factor  of  2  in  the  best-fit  velocity 
with  oxygen  functionalization  implies  that  increase  in  thermal 
transport  is  due  to  an  increase  in  interfacial  bond  strength  by 
roughly  40%. 

We  cannot  preclude  the  possibility  that  the  change  in 
electron  population  from  the  functionalization  could  be 
affecting  the  measured  thermal  boundary  conductance. 
However,  Majumdar  and  Reddy44  showed  that  electron- 
phonon  coupling  could  decrease  the  thermal  boundary 
conductance,  an  opposite  trend  as  what  is  observed  in  our  O- 
SLG  interfaces.  In  addition,  they  concluded  that  the  electron- 
phonon  resistance  only  contributes  to  interfacial  transport 
when  the  phonon-mediated  conductance  is  on  the  order  of 
GW/m2/K,  which  is  over  an  order  of  magnitude  greater  than 
the  values  we  report  here  for  the  thermal  boundary 
conductance  at  our  graphene  interfaces.  Lyeo  and  Cahill45 
experimentally  determined  that  electron  scattering  does  not 
affect  thermal  transport  across  metal/diamond  interfaces.  The 
covalent  nature  of  the  diamond  bond  in  their  study  suggests 
that  the  oxygen  functionalization  of  our  graphene  should  not 
add  any  additional  thermal  mechanisms  via  the  electrons. 
Finally,  Hopkins  et  al.46'47  demonstrated  that  an  increase  in 
thermal  boundary  conductance  due  to  electron  transport  will 
only  occur  during  regimes  of  electron— phonon  nonequili¬ 
brium;  our  TDTR  measurements  are  conducted  during  time 
delays  that  are  orders  of  magnitudes  greater  than  typical 
electron— phonon  nonequilibrium  regimes.  Collectively,  this 
alludes  to  the  fact  that  the  increase  in  thermal  boundary 
conductance  due  to  oxygen  functionalization  is  due  to  changes 
in  the  nature  of  the  chemical  bond  at  the  Al/ graphene  interface, 
and  not  due  to  any  change  in  electronic  contribution. 

The  data  in  Figures  2  and  3  also  give  insight  into  the 
fundamental  phonon  mechanisms  involved  with  thermal 
boundary  conductance  across  Al  and  carbon-based  materials. 
For  example,  without  any  chemical  functionalization,  Al/ 
graphene,  Al/graphite,  and  Al/diamond  all  show  similar  values 
for  thermal  boundary  conductance  despite  the  fact  these  three 
carbon-based  materials  all  have  different  dispersion  relations.  At 
Al/graphene  and  Al/diamond  interfaces,  oxygenation  leads  to  a 
significant  increase  in  thermal  boundary  conductance.  This 
suggests  that  conductance  at  these  interfaces  is  otherwise 
limited  by  interfacial  bond  strength.  For  comparison,  we  plot 
the  ratio  of  our  Al/O-SLG  data  to  our  Al/SLG  data  (Al/O- 
SLG:A1/SLG)  along  with  the  ratio  the  Al/O-diamond  to  Al/ 
diamond  (Al/0-diamond:Al/diamond),  shown  in  Figure  4. 
These  ratios  are  very  similar,  especially  at  elevated  temper¬ 
atures,  which  would  be  indicative  of  a  similar  increase  in  the 
force  constants  at  the  Al— O  bond  compared  to  the  Al-graphene 
or  Al-diamond  bond.  Also  in  this  figure,  we  plot  the  ratio  of  hK 
across  a  Au/(2  nm)Ti/SLG  interface5  (calculated  via  eq  l)  to 
hK  across  an  Au/graphite  interface.14  Presumably,  the  Ti 
wetting  layer  between  the  Au  and  SLG  should  increase  the 
wettability  of  the  Au  to  the  SLG.  However,  the  two 


Figure  4.  Ratio  of  thermal  boundary  conductance  Al/O-SLG  to  Al/ 
SLG  (filled  squares).  For  comparison,  we  also  plot  the  ratio  of  hK 
across  the  Al/O-diamond  interface33  to  the  Al/diamond  interface34 
(filled  diamonds).  We  also  plot  the  ratio  of  hK  across  the  Al/Ti/ 
Graphite  to  the  Al/Graphite  interfacial  conductance  (filled  circles). 
These  ratios  are  very  similar,  especially  at  elevated  temperatures,  which 
would  be  indicative  of  a  similar  increase  in  the  elastic  constants  at  the 
Al— O  bond  compared  to  the  Al-graphene  or  Al-diamond  bond.  We 
plot  the  ratio  of  hK  across  a  Au/ (2  nm)Ti/SLG  interface  to  hK  across 
an  Au/graphite  interface14  (filled  triangles).  The  two  conductances  are 
nearly  identical  over  the  entire  temperature  range  (i.e.,  Au/  (2  nm)Ti/ 
SLG:Au/graphite  ~l).  Where  the  Ti  will  increase  the  adhesion,  Au  is  a 
very  soft  material  with  low  phonon  velocities  and  weak  elastic 
constants.  Therefore,  the  phonon  flux  in  gold  is  so  low  that  regardless 
of  the  strength  of  the  bond  at  the  Au/SLG  or  graphite  interfaces,  the 
phonon  transmissivity  will  always  be  limited  by  the  phonon  flux  in  the 
Au. 


conductances  are  nearly  identical  over  the  entire  temperature 
range  (i.e.,  Au/(2  nm)Ti/SLG Au/graphite  ~l).  On  the 
contrary,  conductances  at  Al:(5  nm)Ti/graphite  interfaces  are 
approximately  a  factor  of  2  higher  than  those  at  Al/graphite 
interfaces  80  to  300  K.  While  Ti  layers  increase  adhesion 
strength  at  metal/graphite  interfaces,  the  increased  adhesion 
only  benefits  the  Al/carbon  interfaces.  This  makes  intuitive 
sense,  as  the  phonon  frequencies  are  much  lower  in  Au  than  Ti, 
and  thus,  are  less  affected  by  a  weakly  bonded  interface.  While 
both  chemical-functionalization  and  the  inclusion  of  an 
adhesion  layer  can  increase  the  conductance  at  Al/carbon 
interfaces,  chemical  functionalization  does  not  require  the 
inclusion  of  additional  materials,  nor  does  it  increase  the 
relative  “thickness”  of  the  interface  (e.g.,  5  nm  of  Ti 
corresponds  roughly  to  the  thickness  of  15  layers  of  graphene). 
Regardless,  the  interplay  between  the  efficiency  of  thermal 
transport  and  the  bond  strength  at  the  metal/SLG  interface  has 
major  implications  for  engineering  and  development  of 
graphene-based  devices  in  which  contact  resistances  plague 
the  device  performance. 

In  summary,  we  have  investigated  the  effects  of  functional¬ 
ization  on  the  thermal  boundary  conductance  across  metal/ 
graphene  interfaces  through  the  introduction  of  chemical 
adsorbates  on  the  SLG  surfaces  in  order  to  increase  the 
density  of  covalent  bonds  bridging  the  metal  and  the  SLG. 
Specifically,  we  metalized  plasma-functionalized  graphene  and 
measured  hK  at  Al/SLG/SiCL  contacts  with  time  domain 
thermoreflectance.  Through  these  measurements,  adsorbates 
on  the  SLG  surfaces  were  shown  to  influence  the  cross  plane 
thermal  conductance.  These  influences  are  attributed  to 
changes  in  bonding  between  the  metal  and  the  SLG.  The 
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thermal  boundary  conductance  was  found  to  be  most 
influenced  by  the  interfacial  chemical  bonding  as  both  the 
phonon  flux  and  the  vibrational  mismatch  between  the 
materials  are  each  subject  to  the  to  the  interfacial  bond 
strength.  In  the  case  of  oxygen  functionalization,  we  show  that 
the  thermal  boundary  conductance  increases,  which  is 
attributed  to  the  increase  in  the  chemical  bond  strength  at 
the  Al/SLG  interface.  The  fact  that  oxygen  functionalization 
leads  to  an  increase  in  thermal  transport  demonstrates  the 
powerful  role  of  chemical  bonding  on  thermal  transport  across 
graphene  interfaces. 
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